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I. INTRODUCTION 



Supersymmetry, the symmetry which links bosonic 
with fcrmionic degrees of freedom, is an intriguing con- 
cept with many applications in quantum field theory and 
statistical physics. It plays a prominent role for open 
challenges in the Standard Model of Particle Physics such 
as the hierarchy problem, and continues to inspire the 
construction of models for new physics. In statistical 
physics, supersymmetry also appears as a technical sym- 
metry in the exploitation of systems which otherwise are 
too difficult to handle. It is thus of great interest to 
further the understanding of interacting supersymmetric 
theories, and to clarify the impact of supersymmetry on 
the phase structure and the critical behavior at lowest 
and highest energies. 

This work is devoted to the supersymmetric extension 
of O(N) symmetric scalar theories in three euclidean di- 
mensions, continuing a line of research initiated in pQ. 
Without supersymmetry, the bosonic theory with a mi- 
croscopic (0 2 ) 3 potential is described by three renormal- 
ized parameters permitting first-order phase transitions 
at strong coupling as well as second order phase transi- 
tions with Ising-type critical behavior [2J. In the limit of 
infinitely many scalars, the analytically solvable spheri- 
cal model also admits an ultraviolet fixed point with bro- 
ken scale invariance, the Bardeen-Moshe-Bander (BMB) 
phenomenon, allowing for a non-trivial continuum limit 
[3H5] ■ With supersymmetry, additional fcrmionic degrees 
of freedom are present and their fluctuations modify the 
quantum effective theory. The O(N) symmetric Wess- 
Zumino model with a microscopic ($ 2 ) 2 superpotential 
is determined by only two renormalized parameters and 
critical and tricritical theories are the same. Its phase 
structure has attracted some attention in the past j6j- 
113] . In the limit of infinitely many superfields, four dif- 
ferent phases have been observed [51 [7], including pe- 
culiar degenerate O(N) symmetric phases with several 
mass scales. Similar to the scalar case, a supersymmet- 
ric version of the BMB fixed point has equally been found 



at a critical coupling where the bosons and fermions be- 
come massive while a Goldstone-boson (dilaton) and a 
Goldstone-fermion (dilatino) are dynamically generated. 
The supersymmetric O(N) model has also been discussed 
in the 1/N expansion [llj . where the authors found a 
non-trivial UV fixed-point and a stable dilaton phase. 
At next-to-leading order the dilaton acquires a mass of 
order 1/N showing that a phase with spontaneously bro- 
ken scale invariance only exists in the limit of infinitely 
many superfields [13] . 

A method of choice in the study of phases and phase 
transitions is Wilson's renormalisation group (RG) |14) . 
It is based on a path-integral representation of the the- 
ory, where the continuous integrating-out of momentum 
modes permits a smooth and controlled interpolation be- 
tween the microscopic and the full quantum effective the- 
ory |15j . Physically- motivated approximations schemes 
together with analytic versions of the RG [TMTB] allow 
for a global analysis and a classification of phase tran- 
sitions and critical exponents even at strong coupling. 
The method has been successfully applied to phase tran- 
sitions in Ising-type universality classes j!5[ I19H21) in- 
cluding high-precision computations of its critical expo- 
nents with increasing levels of sophistication 22 -27] . The 
extension of the functional RG towards supersymmetric 
theories [TJ l2"5r[3"Tj] therefore bears the promise for deeper 
insights into the phases and the critical behavior of su- 
persymmetric O(N) theories. 

This paper is organized as follows: We recall the main 
features of supersymmetric O(N) models including a su- 
persymmetric version of Wilson's RG (Sec. [TT] ) , followed 
by a discussion of its exact analytical solution in the 
large- A limit (Sec. III). We then give a detailed account 



of the phase diagram and phase transitions in the renor- 
malized theory, and examine the appearance of a multi- 
valued effective potential, also in comparison with ear- 
lier findings (Sec. IV). We repeat this exercise with a 
finite short-distance cutoff including a thermodynamical 
derivation of scaling exponents (Sec. [VJ), and examine 
the supersymmetric BMB phenomenon(Sec. VI). At fi- 



nite A, we derive an exact fixed point to leading order 
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in a gradient expansion and evaluate its impact on the 
phase transition, and on the fate of the BMB mechanism 



(Sec. VII ). We close with a brief summary and some con- 
clusions (Sec. VIII). 



II. SUPERSYMMETRIC RG FLOW 

In this section we sketch the features of supersymmet- 
ric O(N) models and recall the supersymmetric renor- 
malization group flow in the local potential approxima- 
tion. For a detailed discussion and derivation see pQ. 

A. Action 

The three-dimensional supersymmetric O(N) models 
are built from TV real superfields 



$i(x, 9) = <j>i(x) + Hi{x) + -99Fi(x) 



(1) 



containing scalar fields </>i, Majorana fermions ipi and 
auxiliary fields Fi as components and a two-component 
anticommuting Majorana spinor 9. The invariant action 



S 



d A x 



~<I>VV<S> + 2W(R) J 



(2) 



wherein we suppress the internal summation index i, con- 
tains the 0(A^)-invariant composite superfield 
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where p = <fi 2 /2. The supercovariant derivatives 

V=-=+i09 and V = -—-i6$ 
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(3) 



(4) 



obey {T>k,T>i} — —21(7^)^/9^. An expansion in compo- 
nent fields yields the off-shell Lagrangian density 

£ g = \ ( 



•-i^ + F 2 ) +W'(p) <t>F 
\W'(p)^~\W"(p){U)(U)- (5) 



By eliminating the auxiliary fields F through their alge- 
braic equation of motion, F = —W'(p) <fi, we obtain the 
on-shell density. The field-dependent fermion mass m^p, 
the bosonic potential V, and the field-dependent Yukawa- 
type coupling Ay all follow from the supcrpotcntial W as 



= W'{p) 

v = p [w'(p)Y 

\y = \W"[p). 



(6) 



All salient features of the classical theory are encoded in 
the functions For a polynomial superpotential the 
scalar field potential always has a minimum at V(0) = 
implying that global supersymmetry is unbroken. 



B. Renormalization group 

Including the effects of quantum and thermal fluctu- 
ations implies that the classical action ^ is modified 
and replaced by a "coarse-grained" or "flowing" effective 
action r k - In the next-to- leading order in the super- 
derivative expansion 



d 3 x 



1 



-<S>Z k WD$ + 2W k 



(7) 



interpolates between the classical action at the high- 
energy cutoff-scale k = A and the full effective action 
at k = 0. The fluctuations above k modify both the 
superpotential, which has turned into a scale-dependent 
superpotential W k , and the kinetic terms, which may ac- 
quire a non-trivial field- and momentum-dependent wave 
function renormalization factor Zk{\$ 2 ,T > T>). 

The RG momentum scale k is introduced on the level 
of the path integral by adding suitable momentum cut- 
offs Rk{q 2 ) to the inverse propagators of the fields. The 
cutoffs regularizes the path integral in the infrared and 
gives rise to a finite flow of the scale dependent effective 
action. Optimized choices for R k are available to ensure 
the stability of the resulting RG equations [16-18 .. The 
scale dependence of the effective action Q is described 
by a functional differential equation 14J 



d t r k 



l. 



-STr T 



.(2) 



dtRk 



(8) 



which emerges as an exact identity from a path inte- 
gral representation. Here, t = lnfc/A denotes the dimen- 

(2) 

sionless RG "time" parameter, 1^ the second functional 
derivative of T k with respect to the fields, and the super- 
trace denotes a momentum integration and a sum over 
all fields, including appropriate minus signs for fermions. 



C. Derivative expansion 

Finally we detail our equations to leading order in a 
super-derivative expansion, the so-called local potential 
approximation (LPA). It amounts to setting the wave 
function factor Z k — 1 throughout, which is a good ap- 
proximation in the large- limit where RG corrections 
to the wave function renormalization of the relevant de- 
grees of freedom, the Goldstone modes, are suppressed 
as 1/N. In scalar O(N) theories, the LPA gives already 
very good results for scaling at the Wilson-Fisher fixed 
point [22 . Here, the LPA does retain the full field- and 
scale-dependence of the superpotential Wh- 
in this work, we introduce the momentum cutoff as 
a supersymmetric invariant F-term of the superfield, by 
adding AS k = \ J d 3 x ^Rk^\§g to the action under the 
path integral, with 



$i? fc (X>X>)$ = -f$r fc (-n)x>x>$. 



(9) 
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The dimensionless function rk(p 2 ) describes the shape 
of the momentum cutoff. The momentum trace is per- 
formed analytically for specific optimized choices for 
[IMS]- Following PUGS], we adopt 



(10) 



The flow in LPA for the superpotential is obtained by 
projecting ^ onto the term linear in the auxiliary field 
F, and this yields 

N fW'\ fW +2pW"\ 

where I(x) = xj (1 + x 2 ). It is understood that W and its 
derivatives are functions of the RG scale k and the fields, 
and we will omit the index k. The first term on the RHS 
is the contribution of the N — 1 Goldstone modes and the 
last term is the contribution of the radial mode. Note 
that the RHS of the flow vanishes for W' = 0, and for 
1/| W'\ — > 0, corresponding to the classical limit where 
the couplings and the potential §6§ are independent of 
the RG scale. 



To achieve the simple form (111 we have rescaled the 



fields and the superpotential as 



N N 
p, W-> — W. 



8tt 



8tt 2 



(12) 



Note that W is invariant under the rescaling which ab- 
sorbs the redundant overall factor l/(87r 2 ), originating 
from the momentum integration, into the field and the 
superpotential. The additional rescaling with N also re- 
moves the leading ./V-dependence from the RG equation 



(111. In these conventions, and with given initial condi- 
tion Wk=K (p) the RG flow determines the superpotential 
in the infrared limit k — > 0. 

To study the critical behavior we introduce a dimen- 
sionless field variable p, a dimensionless superpotential w 
and a dimensionless scalar potential v as 



P 
k ' 



w(p) 



W{p) 



v{p) 



W'(p) 



.(13) 



In terms of ( 13 1 the flow equation ( 11 ) reads 



d t w+2w-pw' = -(l-^)l(w')-^I{w'+2pw") . (14) 
For completeness we add the flow equation for w' = u, 



1 



d t u + u-pu' = -(1- — )u'I'(u) 



N 



(3vf + 2pu")l'(u + 2pu') , (15) 



and similarly for higher derivatives of the superpotential. 



III. EFFECTIVE POTENTIAL 

In this section, we discuss the explicit and exact solu- 
tion for the effective potential in the limit 1/N — > 0, and 
derive the main equations which govern the symmetry 
breaking in this model. 



A. RG flow and boundary condition 



In the large- N limit, the flow equation ( 15 1 for u = w 
simplifies considerably and is given by 



d t u 



1 



pu 



(1 



(16) 



The terms on the LHS encode the canonical scaling of 
the superpotential and the fields and the RHS encode 
the effects due to fluctuations. The integration of ( 16 1 



with respect to the logarithmic RG scale t = fnfe/A gives 

F{u) = G(ue*) (17) 



P-1 

u 



with 



F(u) 



1 



+ 2 arctan(u) 



(18) 



The function G(x) is determined by the initial conditions 
for u{p) imposed at some reference scale k = A. We use 
throughout the boundary condition 



k = A 




t(p-k) 

T {p — K A) 



(19) 



where r denotes the quartic superfield coupling at the 
cutoff. We recall that it is an exactly marginal coupling, 
i.e. that dtT = 0. If the UV parameter n is positive, n A 
is interpreted as VEV for the scalar field at k = A. 

Following CQ, the fixed point solutions are 
parametrized in terms of the parameter 

c=1/t. (20) 

Then the function G(x) is given by 

G{x) = c - F(x) + 

in terms of the initial parameters. For initial conditions 
different from ( 19 1 the function is modified accordingly. 



K- 1 



(21) 



B. Factorization 



Using the initial condition ( 19 1, the analytical solution 
(17) takes the form 



p-po(t) =cu + H(u)-H(ue t )e~ t 
p (t) =l + SKe~ t , t = hik/A, 
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where the non-negative function 



identified: 



H{u) =uF{u) 



1 



2u arctan u 



(23) 



encodes the RG modifications due to fluctuations pQ. 
The parameter 5k — k — 1 measures the deviation of 
the VEV at the initial scale po(t = 0) = K from its crit- 
ical value k c1 - = 1. For any positive deviation we have 
Po(t) — > oo in the infrared limit corresponding to a fi- 
nite VEV of the scalar field. Since the potential V in 
^ shows a second minimum at p = 0, the global O(N) 
symmetry is (not) spontaneously broken if the finite (van- 
ishing) VEV is taken. Conversely, for a negative 5k we 
have po(t) < in the infrared limit such that the global 
minimum of the effective potential is achieved for vanish- 
ing p. This leaves the global O(N) symmetry intact. The 
case 5k = then corresponds to the boundary between 
the symmetric and broken phases. 

From ( 22 ) we conclude that the IR repulsive mode 



associated with po(t) is solely controlled by the initial 
VEV, independently of the coupling strength r. This 
has been seen previously in purely scalar theories in the 
large- TV limit [2D]. All the remaining couplings included 
in the potential are either exactly marginal or IR attrac- 
tive. Their flow is encoded in the term H(ue)e in the 
first equation of (22 1. This factorization of the solution 
is a consequence of the large-TV limit, and allows for a 
straightforward analysis of the entire phase structure of 
the model. The global form of solutions u(p, t) is mainly 
determined by the coupling r = 1/c and the function H, 
with po only entering through a shift of the p-axis. 

The non-negative function H appearing in the implicit 
solution ( [22] ) will be of importance below. Expanding H 
in powers or 1/u leads to 



H = n\u\-1 



1 

3^2 



1 



Conversely for small u we find the expansion 

H = 3u 2 --u i + -u 6 + 0(u s ). 
3 5 



(24) 



(25) 



The solution (22 1 is invariant under (c, u) -H- (— c, — u) 
since H(u) is an even function. Furthermore, the scalar 
field potential only depends on u 2 and we may restrict 
the discussion to c > 0. 



C. Fixed points 

We briefly recall the main results from [I]. The fixed 
point solutions follow from (171 by setting G(ue*) to a 
constant c, 



p = l + H(u*) 



(26) 



The constant c is related to the marginal quartic super- 
field coupling t = u'(p = Po(t)) as c = 1/r. Five charac- 
teristic values ci < cl < cp < cm < cq for |c| have been 



ci 
cl 
cp 

cm 
cg 





Ktt + 3) 
n 

§V3 



(27) 



The extreme values cj and cq correspond to the 'would- 
be' Wilson-Fisher and the Gaussian fixed point solution 
p{u^) 1 respectively. The solutions exist and extend over 
all physical field space p > in the weak coupling regime 
c p < I cj < cg- For \c\ > cm, fixed point solutions are 
monotonous functions of and extend over the entire 
real axis. In the intermediate coupling regime cl < |c| < 
Cp, fixed point solutions exists both with and without a 
node at p = 1. Finally, in the strong coupling regime 
I cj < cl, the solutions do not extend over all fields p > 0. 
Numerically, the ranges 



cm - cp 
cp 



0.011, 



cp - c L 
cp 



0.023 



(28) 



are very small. The fixed points are non-Gaussian except 
for |c| = cq, yet they displays Gaussian scaling for all 
physical fixed points except for jcj = cp or cj. 



D. Non-analyticities 

Finally, we discuss the appearance of non-analytic be- 
havior in the integrated flows at intermediate and strong 
coupling. This discussion completes the general descrip- 
tion of fixed point solutions in [T] and will be of help to 
understand the RG flows away from critical points in the 
next section. 

By construction, the basic flow equation ([8| is well- 
defined (finite, no poles). Furthermore, the RHS of the 
supersymmetric flow ( 16 1 is bounded, provided that the 
superpotential remains real. Incidentally, this is in con- 
trast to the standard purely bosonic flows, which po- 



tentially may grow large in a phase with spontaneous 
symmetry breaking. Despite their boundedness, the su- 
persymmetric fixed point solutions display Landau-type 
poles at strong coupling due to non-analyticities, such as 
cusps, of the integrated RG flow. This can be appreci- 
ated as follows: consider the field-dependent dimension- 



less mass term u\p). From the fixed point solution (26) 
we conclude that it diverges provided that 



dp I 
du ' 



= c + H'(u s ) = 0. 



(29) 



This condition determines the singular value u s and from 
(26) we obtain the value of the singular field, 



l + H(u s )-u s H'(u s 



1 



(l + u2)' 



(30) 
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FIG. 1: The non-monotonic odd function H'(u). 



The function H'(u) is odd and bounded by H'(u c ) = 
±cm- Asymptotically, we have \H'(u — > ±oo)| = cp < 
cm, see Fig. [I] Hence, with decreasing |c| a divergence 
for u! is first encountered for |c| = cm- Performing an 
expansion of ( |26[ ) up to the first non-trivial order, we find 
that 



p - Pc = ^-ff"'(tt c )(u» - u c f 



(31) 



up to subleading corrections. In the expansion we used 
(29) and that H" vanishes at u c . We note that (31) is 



continuous across (u*,p) = (u c ,p c ). Therefore, the non- 
analyticity in the solution can be written as 



T sgn(p - p c 



P- Pc 



\H"'(u c 



1/3 



(32) 



where the signs refer to c = Tcm, leading to a non- 
perturbative Landau pole in u' t , 



P~ Pc 



,2/3 



(33) 



At a fixed point solution, the Landau pole remains invis- 
ible, because it is achieved at the negative p c — —1/8. 
Increasing the coupling by lowering |c| below cm, the ex- 
pansion in the vicinity of dp/du — becomes 



p - Ps = )^H"(u s )(u* - u s ) 2 



(34) 



up to subleading terms, where u s is determined through 
(29). In this regime, H"{u s ) is non-zero throughout. In 



the parameter range cp < \c\ < cm we find two solutions 
for u s with \u s i\ < \u c \ < \u S 2\ and H"(u s i) < < 
H"{u S 2)- Effectively, the solution for the superpotential 
becomes multi-valued in a limited region of field space. 
For |c| < cp we find one solution for u s with H"(u s ) > 0. 



In contrast to (33), the non-analyticity has turned into a 
square root, 



1 



±2\H"(u s )\ 1/2 (p-p 



,1/2 



(35) 



The non-analyticity (|35j) is stronger than (|33j) and the 
solution ( 34 ) cannot be continued continuously beyond 



the point (u*,p) = (u s ,p s ). For |c| < cl, we have that 
Ps(c) > and the pole appears in the physical regime. 
In contrast, the solutions extend over all fields provided 
that p s < which is the case for |c| > cl- 

It is interesting to note that non-analyticities, such as 
cusps, have been detected previously in the context of the 
random field Ising model, where disorder is technically in- 
troduced with the help of Parisi-Sourlas supersymmctry. 
Using functional renormalization, it has been argued that 
a cusp behavior at finite "Larkin scales" k = k^ > is at 
the origin for the spontaneous breaking of Parisi-Sourlas 
supersymmetry [4*0H4"2"] . 

At this point it should be mentioned that the superpo- 
tential W' shows another non-analytic behavior: It is not 
differentiable at its node po in the exact IR limit for ar- 
bitrary couplings c > 0. This issue is discussed in detail 
in Sec. IVB21 and |VD] below. 



IV. RENORMALIZED FIELD THEORY 

In this section, we discuss the spontaneous breaking of 
symmetry and the phase structure of the model in the 
limit where the UV scale A is removed. 



Renormalization 



The solution (22) is valid for all k and A, and we may 
take the 'continuum limit' 1/A — > 0. The term contain- 
ing the explicit i-dependence drops out in the continuum 
limit, in consequence of the limit k/A — > for fixed and 
finite k and (25). The remaining scale-dependence solely 



reduces to the implicit scale-dependence of po(k) in 



Po(k) 
Po(k) 



cu + H(u) 
1 + Po/k . 



(36) 



The dimensional parameter po has taken over the role of 
5k A in ( 22 ) . In the above, the VEV (or the mass term, 
respectively) is the only quantity which is non-trivially 
renormalized in the continuum limit by requiring that 



po = lim (Sk(A) A) < oo . 



(37) 



Consequently, the canonical dimension of fields remain 
unchanged (no anomalous dimension). The continuum 
limit maps the original set of free parameters (t, k, A) to 
the parameters (r, po). Note that all couplings of the su- 
perfield derivative - the marginal coupling c and the IR 
attractive higher-order couplings u^(po) - have settled 
on their fixed point values. The only 'coupling' which 
has not settled on a fixed point is the UV attractive di- 
mensionless VEV po. With this perspective, p and the 
non-renormalized parameter c should be viewed as a free 
parameters of the model, fixed by the microscopic pa- 
rameters of the theory. In terms of the dimensional fields 



G 



p = pk and superfield derivative W'(p) — 
integrated RG flow becomes 

p-po(k) =cW' + kH(W'/k) 
p (k) =k + p . 



u(p) k, the 



(38) 



We note that po also has the interpretation of the physical 
VEV in the infrared limit of the theory, provided it is 
positive. Below, we find it is useful to switch between 



the representations (36) and (38 1 



B. Characteristic energy 



The RG flow (36), (38) carries a characteristic energy 
scale E, meaning that the theory changes its qualitative 
behavior depending on whether fluctuations have an en- 
ergy larger or smaller than E. The energy scale is set by 
the UV renormalization of the model (37) and given by 



E=\po\ 



(39) 



For k 3> E, the dimensionful VEV scales proportional 
to k, and the dimensionless parameter po becomes a con- 
stant. This corresponds to a fixed point. All other dimen- 
sionless couplings equally have stopped to evolve with 
RG scale and thus the entire solution approaches a high- 
energy (UV) fixed point. This fixed point would persist 
for all k provided that E = 0. It then has also the in- 
terpretation of an IR fixed point. This regime is most 
conveniently described using (36). For E > 0, and with 



decreasing k, deviations from the fixed point become visi- 
ble once k reaches E. Here the VEV displays a cross-over 
from linear scaling p{k) oc k for k 3> E to the constant 
value po for k <C E. In full analogy, the dimensionless 
VEV displays a cross-over from a constant value to scal- 
ing inversely proportional to the RG scale. In addition, 
the running of all dimensionful couplings in the poten- 
tial is switched on once k ~ E and below. This regime 



is conveniently described using ( 38 ) which governs the 



remaining RG running through its RHS. 



C. Gap equations 

We first discuss the phase structure based on the in- 
tegrated RG equations in the IR limit k = 0, see Fig. [2] 
This allows for a direct comparison with earlier results 
based on gap equations and Schwinger-Dyson equations 
[H [7] . In the infrared limit we may use ( 24 ) in ( 38 ) and 
obtain 



po = cW' + c P \W'\ 



(40) 



Since the potential shows a local minimum at vanishing 
field, the squared particle masses are given by 



fj, 



v"m <f> = = w' 2 (p)\- p=0 . 



(41) 



cp 



SYM 


SSB 


m 


M M p 


m M 





FIG. 2: Schematic phase diagram of the supersymmetric 
model based on the gap equation ( |42[ ) in the infinite cutoff 
limit. Results agree with earlier findings in [51 [7J. 



Thus, (40) becomes a gap equation for the mass param- 
eter ft = W'(p 



0), 



Po 



-cp - c P \p\ 



(42) 



The significance of (42) is as follows. For fixed p and c 
it yields the possible infrared solutions for the masses at 
vanishing field. Without loss of generality we restrict the 
discussion to c > 0. For non- vanishing po we find two 
solutions 



m = p 
M= -p 



po 



Cp + c 

Po 



cp 



> 

> 0. 



(43) 



In the symmetric regime with negative po the mass m 
is always present and the second mass M is available as 
long as c < cp. In the SSB regime with positive po there 
are two degenerate ground states: As expected, we find 
a non-symmetric ground state with a radial mass M t 



see sections 



pi 



V D and VE However, for c > cp, the gap 
equations show an additional symmetric ground state, 
characterized by the mass M. Note that changing the 
sign of c leads to equivalent results under the following 
replacements 



(c, m, M, M p ) o (-c, M, m, -M p ) . 



(44) 



At the phase transition, i.e. for po = 0, the gap equation 
( |42| ) states that either c = tt with the mass M > un- 
determined, or c = — 7r and the mass m undetermined. 
These findings agree with the earlier ones from [C][7|. The 
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this paper po r = 1/c 

Bardeen et.al. [BJ -A-k 2 ^ 1 (4tt 2 )" 1 A 

Moshe and Zinn- Justin [Jj — 47r 2 (p — ^ c )it _1 (4tt 2 )^ 1 u 



TABLE I: "Translation guide" between the conventions used 
in [BJ, [7j, and this paper. 



sole difference is that the value for the critical coupling, 
cp, depends on the regularization. The precise link to 
the conventions used in [7] is given in Tab. pi 



D. RG phase diagram 

Next we discuss the phase diagram implied by the in- 
tegrated RG equations for all scales k, and compare with 
the results based on the k = limit. 
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1. Graphical representation 

We begin with a useful graphical representation of the 
rcnormalized RG trajectories (36). For vanishing p~o, we 



note that the trajectories (36) reduce to the fixed point 
solutions u*(p) analyzed in [I]. The only difference with 
the fixed point solutions is related to a shift of the argu- 
ment, 

u(p)=u4X), X = p+l- Po (k) = p- P ^ (45) 

in terms of the fixed point solutions. 

The structure of the solutions and their dependence on 
the constant c is shown in Fig. [3j Once the free param- 
eters are fixed, the RG evolution of a particular solution 
stays on a curve with constant c, indicated by the curves 
given in the Figure. Rotating counter-clockwise around 
(X,u#) = (1,0) from the horizontal cg-line to the cj- 
curve (from the c/-curve to the cc-line) covers all curves 
with positive (negative) c. Both sets connect through the 
point (1,0). We recall that (c, it*) <-» (— c, — u*) describe 
equivalent physics. 



Using (45) and (36), we conclude that for u(p) to cover 
all physical fields p £ [0, oo], we need that 



X g [-po/k,oo] 



(46) 



The curves it* (X) in Fig. [3] define monotonous (and in- 
vertible) functions provided that X > 1. A unique classi- 
fication of curves is then achieved by choosing a value for 
u„ on a line of constant X > 1, together with fixing po. 
Interestingly, two different values for may correspond 
to one and the same parameter c. Below, we mostly stick 
to the classification in terms of c, and we will highlight 
situations where this is no longer sufficient. 



FIG. 3: Graphical representation of the solutions u*(X) of 
(36), where X — p — po/k. The shaded areas are separated 
by thick lines at \c\ = cj, c L , c P , c M and cg- 



2. Symmetric regime 

The symmetric phase is characterized by a finite and 
negative p~o and for large scales X reduces to p. A re- 
striction on the coupling parameter c is imposed if we 
require that the solution u should exist for all p. For 
weak coupling, 



Po < , \c\ > c P 



(47) 



all it* are single- valued for non-negative arguments such 
that the u*(X) stay well-defined for all scales, see Fig. [4] 
For intermediate coupling 



po < , c P >\c\> c L 



(48) 



the theory admits two distinct effective potentials, and 
two scalar mass parameters. They are related to trajec- 
tories which either run through a node, or not, depending 
on whether u(0) for k E is larger or smaller than 1, 
see Fig. |4j The theory is then characterized by the cou- 
pling and the scalar mass at vanishing field. This peculiar 
structure has been found previously and we discuss it in 
more detail below. 



3. Symmetry broken regime 

Spontaneous symmetry breaking is possible for positive 
Po- This requires that u*(X) has to be well-defined for 
all real X. In view of the analytical solution in Fig. |3j 
this limits the achievable couplings to 



po > 0, | cj > c P . 



(49) 
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FIG. 4: RG trajectories in the O(N) symmetric phase: at 
weak coupling, trajectories either show a non-vanishing VEV 
for large scales (SYM, yellow shading), or a vanishing VEV 
for all scales (SYM, green shading). At strong coupling tra- 
jectories terminate at Landau poles. 



FIG. 5: RG trajectories u*{X) according to ( |36[ ) in the re- 
gions with spontaneous breaking of the O(N) symmetry in 
the parametrization (451. Couplings are either finite for all k 



(SSB, blue shading), or run into a singularity (Landau, red 
shading). Some trajectories cannot be continued beyond the 
Landau pole (magenta shading). The SSB phase cannot be 
defined for strong coupling (white area). 



Smaller |c| do not lead to a well-defined physical theory 
in the IR. For 



Po > 0, |c| > c M 



(50) 



the function is one-to-one and the theory described 



by u(p) in (36) remains well-defined even in the IR limit. 
The theory is then characterized by two mass scales. The 
first one is given by the scalar mass at vanishing field 
corresponding to an O(N) symmetric phase, whereas the 
second mass scale is given by the radial mass at p — po 
allowing for SSB. 



4- Strong coupling and Landau regime 

It remains to discuss the strong coupling and Landau 
regimes in Figs. [4] and [5j We begin with trajectories in 
the SYM regime, with 



Po < 0, 



c < c L . 



(51) 



We take a 'bottom-up' view according to which the cou- 
plings evolve from the infrared towards higher scales, 
parametrizing the effective potential in terms of local 
couplings in an expansion about vanishing field. Tra- 
jectories with (51) emanate from the upper /lower- right 
corner in Fig. |4j for k « and increasing k corresponds 
to decreasing X. With increasing fc, the running mass 
term and the fermion-boson coupling at vanishing field 



u' {p = 0) = u'^—po/k) diverge at k — k^, and the rcnor- 
malized RG flow comes to a halt: the solutions p6| ) can- 
not be continued beyond these points, because X cannot 
decrease any further along the integral curve u*(X). In- 
terestingly, the potential is double- valued for k < kr, with 
two different trajectories terminating at the same Landau 
pole. Using (34) together with (45), the non-analyticity 
in u reads 



P-PM = \H"(u s )(u{p)-u s f 
and the Landau poles are located at 



p s (k) = p 3 - 1 + Po(k) = p s + 



Po 

k 



(52) 



(53) 



From the fixed point solution we know that p s < 1 and 
therefore p s (k) < po{k) for all k. In the IR limit, this 
implies that p s (k) — > po(k) from below. Here, the values 
for p s are fixed by the coupling strength c via (30) and 
is positive in the regime ( [51] ). From (53) it follows that 
ki = —p /p s is positive, see Fig.[6j We conclude that the 
parameters (51) allow for a supersymmetric model with 
linearly realized O(N) symmetry up to scales k = k^. 

Next we discuss the SSB regime starting with interme- 
diate couplings 



po > 0, c P < \c\ < c M ■ 



(54) 



Here all curves of constant c contain two Landau points 
with |w s i| < \u c \ < \u s2 \ and H"(u s \) < < H"(u s2 ) 
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FIG. 6: Location of the UV Landau pole for trajectories in 
the SYM phase at strong coupling with (51 1 (see text). 



(see Sec. HID). Each of them is described by (52) with 
(53) and parameters > p S 2{c) > p s i(c)- The singular- 
ity at (p s i,u s i) corresponds to an IR Landau pole ('top- 
down'), whereas the one at (p S 2,u S 2) corresponds to an 
UV Landau pole ('bottom-up'). In the infrared limit, the 
domain where u is multi-valued, collapses to a point with 
Ap = k (p s i — p S 2) — > 0. The location of both disconti- 
nuities approach the VEV p s (k) — > po from below, and 
the discontinuity in the superpotential derivative 

AW' = W'(p s i) - W'(p s2 ) = k(u sl - u s2 ) (55) 

then also becomes arbitrarily small. Interestingly, the 
UV and IR Landau poles become degenerate on the in- 
tegral curve for |c| = cm where p s i — p s2 = —1/8. The 
non-analyticity evolves with 



u(p) 



Tsgn(p- p s (k)) 



P~ Ps(k) 



-H" 



(u s ) 



1/3 



(56) 



together with (53). In this case, the quartic scalar self- 
coupling u'{p = 0) still diverges at the Landau pole, but 
the renormalizcd RG flow continues non-perturbatively 
rendering u'(0) again finite. The non-analyticity (56) 
first appears for vanishing field at the scale k^ = —po/Ps 
and evolves up to the VEV po in the IR limit. 

Next we consider the SSB regime at strong coupling, 



Po > 0, 



c < c L . 



(57) 



The model has a radial mass proportional to the VEV. 
Curves of constant c in Fig. [5] have a Landau pole with 
( 52|53 ) and parameter p s > 0. The integral curves have 
no continuation beyond the pole, which occurs within 
the physical regime for all k. In particular, the effective 
potential is not defined for the entire inner part p < po in 
the IR limit and a scalar mass IF' (0) cannot be defined. 

Finally we consider trajectories in the SSB regime, 
with 
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cp 
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m 
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M M p 


Landau 
M M p 
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strong coupling 
Mp (p > po) 
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m M (k < k L ) 





FIG. 7: Schematic phase diagram based on the RG in the in- 
finite cutoff limit. The scale fcz, is given in Fig. [6] The tighter 
constraints as opposed to Fig. [2] arise from the inspection of 
the full effective potential at all scales k. The parameter range 
between cm, cp and cl is very narrow (281. 



cover all positive 
. In a 'top-down' 



Here, in contrast to (57), solutions 
values for X even for large k; see i 
perspective (with decreasing k) trajectories in the regime 
( [58] ) emanate at X w and continue towards smaller X. 
Again, all trajectories reach a Landau pole for the quartic 
(and higher) superficld coupling at vanishing field, given 
by ( p52| and ( [53] ) with the parameter p s (c) taking nega- 
tive values. The Landau scale reads k^ — —p /p s > 0, 
and the effective potential does not exist for fields below 
p s (k) = k(p s — 1) + po(k) < Pa{k). As in (57), the the- 
ory still has a radial scalar mass set by the VEV and the 
quartic coupling, because the one-sided derivative ^%-\p 
with p> po can be taken for fields larger than the VEV. 
In turn, a scalar mass at vanishing field cannot be de- 
fined. Therefore we conclude that the renormalized RG 
flow cannot be continued towards the infrared for scales 
below the Landau scale k < k^ for parameters ( 58 1 . 



po > 0, c L < \c\ < c P . 



(58) 



E. Discussion 

Our results are summarized in Fig. [7] and should be 
compared with Fig. [2] The phase diagram is given in 
dependence on the coupling parameter c and the scale 
parameter po- 

In the SYM regime, the theory has a weakly coupled 
phase with a scalar mass m where both the O(N) sym- 
metry and supersymmetry are preserved ( |47[ ). With in- 
creasing coupling parameter r, the theory admits two 
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O(N) symmetric phases with two mass scales m and M 
( 48 1 . This regime has a very narrow width in parameter 



space, see (28), which is sensitive to the underlying regu- 
larization. For strong coupling (51), the theory displays 



two mass scales m and M. However, it is also plagued by 
Landau-type singularities which admit no solution for the 
superpotential at scales above the Landau scale fcj,. This 
is not visible from an evaluation of the IR gap equations 
alone, see Fig. [2] for comparison. 

In the SSB regime, the theory has a weakly coupled 
phase (50) where the O(N) symmetry could be sponta- 



neously broken and the effective potential for the scalar 
has two degenerate minima corresponding to two mass 
scales M and M p . The first mass scale is associated 
with an O(N) symmetric phase, whereas the second mass 
scale emerges from a finite VEV allowing for SSB. Fur- 
thermore, global supersymmetry remains intact. With 
increasing coupling parameter r, the theory enters a nar- 
row parameter range where RG trajectories would run 
through a series of Landau poles at intermediate ener- 
gies (54). Here, the discontinuity in field space and in 



the superpotential derivative shrinks to zero in the IR 
limit, the details of which are sensitive to the underly- 
ing regularization. For even larger couplings |c| < cp 
( 57 ) and ( 58 ) , the theory is so strongly coupled that RG 



trajectories terminate at Landau poles in the physical 
regime. The effective potential does not exist for fields 
below the non-trivial VEV p < po in the IR limit. Still, 
the potential does admit a radial mass M p . 

Unbroken global supersymmetry requires a ground 
state with vanishing energy, and an elsewise positive di- 
mensionful effective potential for all fields and all RG 
scales. Strictly speaking, the non-existence of an effec- 
tive potential for small fields means that we cannot de- 
cide, based on the potential alone, whether supersymme- 
try is spontaneously broken at strong coupling, or not. 
However, the occurrence of a Landau scale k^ makes it 
conceivable that supersymmetry may be spontaneously 
broken in the strongly coupled regime. This interpreta- 
tion would be consistent with the picture for the sponta- 
neous breaking of Parisi-Sourlas supersymmetry in dis- 
ordered Ising models @0] , which is triggered by cusp-like 
non-analyticities of the RG flow at a finite " Larkin scale" 
k,L. At strong coupling, these limitations of the full effec- 
tive potential and the occurrence of Landau poles are not 
directly visible from the infrared limit only, see Figs. [2] 
and [7} It is a virtue of the fully integrated RG flow for 
all scales k that the structure of the effective potential at 
strong coupling has become transparent. 



V. EFFECTIVE FIELD THEORY 

In this section we discuss the integrated RG flow from 
an effective theory perspective. We assume that the UV 
scale A is finite, and that the boundary condition at k = 
A has been achieved by integrating-out the fluctuations 
with momenta above A. The RG equations then detail 
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FIG. 8: The graphical solution of the gap equation (611 in 
the symmetric phases with po < 0. For positive p there exist 
one, two or three solutions. 



the remaining low-energy flow of couplings for all scales 
k < A. In terms of dimensional quantities, the solution 



(22) reads 



p-p (k) =cW' + H 
p a (k) = k + po . 



W 

IT 



k - H 



W 

X 



A 



(59) 



The parameter po is given by po = A(k — 1) in terms 
of the microscopic (UV) parameters. Our motivation for 
studying (59) is twofold. Firstly, we want to further clar- 



ify the origin of the "peculiar" phases discussed in the 
previous section. Second, we want to evaluate the effect 
of changes in the boundary condition and higher-order 
couplings on the phase structure and critical phenomena 



A. Gap equations 

We begin with the IR limit of the integrated RG flow. 
The corresponding gap equations for the scalar masses at 
vanishing field W'(0) = [i A, i.e. in the O(N) symmetric 
phases, are given in terms of the dimensionless parameter 
p by 



pa 

'A 



cp + c P \p\ - H(p). 



(60) 



where we used expansion (24) for H. For po ^ we find 



two possible branches of solutions with 



H(p) =(c + c p )p+^, ( M >0) 



H{p) = {c-c p )p + 



A 
po 
A' 



(p < 0). 



(61) 



We consider c > since changing the sign amounts to 
interchanging p -H- — p in ( |61[ ). The main difference with 
( 42 ) in the infinite cutoff limit is the appearance of the 



term H(p). 
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FIG. 9: Phases of the supersymmetric model according to the 
gap equations at finite UV scale. The SYM phase displays 
either a single mass scale m, or several ones. The SSB regime 
displays two scalar mass scales M and M p for all couplings. 
Note that the mass scale M represents an O(N) symmetric 
state within the regime where we would normally only expect 
SSB to occur. The phase transition between the SYM phase 
and the SSB phase is continuous with Gaussian exponents. 
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FIG. 10: Phases of the supersymmetric model according to 
the RG equations at finite UV scale. The SSB regime is quite 
similar to the result at 1/A = 0, see Fig[7] The SYM phase is 
substantially larger (see text). The phase transition between 
the SYM phase and the SSB phase is continuous with Gaus- 
sian exponents. Note that there exists a very tiny Landau 
phase for couplings cp < jcj < cm in the SSB regime (red 
line). 



In the SYM regime with negative po we find one, two, 
or three solutions to (61) with m — A/i > 0, and none, 
one, or two solutions M = —A/i > 0, see Fig. [8j Three 
solutions for positive A/i can only exist if the slope c + cp 
is inbetween cp and cm, cf. Fig. [TJ 

For most parts of the parameter space we only have a 
single scalar mass m, similar to the weak coupling phase 
of the renormalized theory. For small po/A and strong 
coupling, a "triangle" opens up allowing for two addi- 
tional mass scales of the type M. The borderline c(po/A) 
is found analytically, starting at the point (c, po/A) = 
(c P , 0) and ending at (c, po/A) « (0,-1.077), see Fig. [9] 
Furthermore, we find two more masses of the type m in a 
tiny "spike" -like region at very strong coupling, bordered 
by the curves connecting (c, po/A) — (cm — cp, — 9/8) « 
(0.035,-1.125) with (c,p /A) = (0,-1.077) and (0,-1) 
as indicated in the same Figure. The bordering lines 
c(po/A) are known analytically. In total, we either have 
a single mass m, or three masses m + 2M or 3m, or five 
different mass scales of the type 3m + 2M in the region 
where the triangle and the spike overlap. Some of the 
masses are parametrically large in the strong coupling do- 
main. We believe that these masses in the very strongly 
coupled domain are an artifact of the regularisation and 
should not be trusted. 

In the regime po > allowing for SSB, a unique scalar 
mass solution M to (61 1 is achieved from the branch with 



negative p, for all couplings. In addition, the theory 
shows the expected radial mass M p . However, we em- 
phasize that some of the solutions found here, in partic- 
ular those at strong coupling, have parametrically large 
masses suggesting that these may be spurious. 



B. RG phase diagram 

Next we turn to the phase diagram of the integrated 
RG flow at finite A for all scales k. 



1. Symmetric regime 



The phase diagram corresponding to (59) is given in 
Fig. |10| where the axes denote the (inverse) quartic su- 
perfield coupling 1 jr and the parameter po in units of the 
initial scale A. For po < the theory is in the symmetric 
phase, provided that the coupling is small enough. There 
is also a strong coupling regime where the RG flow de- 
velops a Landau pole and the effective potential becomes 
multi-valued in the physical regime p > 0. The boundary 
between the two regimes is marked by a curve c cr (po /A) . 
The latter is determined as follows: In the IR limit, the 
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FIG. 11: Possible values W' L /A as a function of the inverse 
superfield coupling coupling c = 1/r associated with a Landau 
pole in the IR limit, i.e. with dp/dW' = 0. 



solution (59) reads 



p- p = cW' +n\W'\- H 



w 



A 



(62) 



and shows a Landau pole, if dp/dW' vanishes. Using ( 62 ) 
in the condition dp/dW' = at p = yields 



Po 
A 



n H , (w L 



A 



A 



(63) 



where W£ is equal to W'(0) when the Landau pole en- 
ters the physical region at p = 0. The real roots of this 
polynomial equation are 



Wl 
A 



± 



y/9 + 8p /A - (3 + 2p /A) 
2 (1 + Po/A) 



1/2 



(64) 



where the plus (minus) sign belongs to the critical cou- 
pling characterizing a Landau pole at p = in the posi- 
tive (negative) half-plane of W. Inserting this into (62 1, 
evaluated at p = 0, yields the critical couplings 



1 

w 



-Po 



H 



W 

X 



A — c P \W' 



(65) 



iri 



as a function of the VEV p~o. In general, we find that 
the occurrence of Landau poles is only possible in the 
parameter range 1 p~o/A € (—1.125,0) and c £ (0,cp), 



i.e. the strong coupling regime. Besides ambiguities with 
W' L < for couplings c < cp, we also find Landau poles 
with W' L > in the very narrow strong coupling regime 



with c < (cm — cp) ps 0.035, see Fig. 11 



Hence, we interpret the different regimes of the sym- 
metric phase as follows (see Fig.JlO]): We observe Landau 
poles in the physical regime with W' L < 0, if the super- 
field coupling r is larger than c~_, i.e. c < c cr _. The 
corresponding borderline starts at the point (c, p~o /A) — 
(cp,0) and ends at (0,-1.077), similar to borderline re- 
sulting from the gap-equation analysis. Furthermore, for 
very strong couplings c < cm - cp < 1 we observe ambi- 
guities with W'j > in the physical regime (dark shaded 
area in Fig. |10[ ). However, this area is bounded by c cr+ 
from below, where c cr + starts at (cm — cp, —1.125) and 
ends at (0,-1.077). 

Interestingly, the available domain of couplings is sub- 
stantially larger than in Fig. [7[ The reason for that is 
quite intuitive, since decreasing the VEV po comes along 
with a shift of the solution W to the left and thus the 
Landau pole may enter the unphysical regime p < 0. In 
addition, the equations do not admit a second mass M, 
unlike the case for 1/A = 0. We emphasize that the 
RG study of the phase diagram also allows for a simple 
descriptive explanation of the occurrence of the various 
masses as shown in Fig. [9j The two additional masses 
M, observed in the strong coupling domain (see big tri- 
angle, Fig. [9]) result from an ambiguity of the solution 
W' in the negative half-plane. The borderline connect- 
ing (cp, 0) and (0, —1.077) in Fig. [^represents the special 
solution with c = c cr _ showing a Landau pole in the IR 
exactly at p — and this corresponds to an additional 
infinitely large mass M. Similarly, the two additional 
masses of type m in the spike-like strong coupling region 
result from ambiguities of the solution for positive W. 



2. Symmetry broken regime 

For po > 0, the theory is in a phase featuring spon- 
taneous O(N) symmetry breaking. For sufficiently weak 
coupling with |c| > cm, the theory displays a well-defined 
low-energy regime with two mass scales M and M„. The 
first one is associated to the curvature at vanishing field 
and thus represents an O(N) symmetric phase, whereas 
the second mass is given by the curvature at the non- 
vanishing VEV po arL d implies SSB. 

In the very narrow coupling-regime cp < \c\ < cm 
there occur IR Landau poles at (p s (k),u s (k)) with 



p s {k) = 1 + cu s + H(u s ) - H (u a e*) e~* + 
= c + H' (u s )-H' (u s e 4 ) 



Po 



(66) 



1 Note that we allow for negative n = po/ 'A + 1, i.e. classical po- 
tentials with a single minimum at p = (symmetric phase). 



within the physical regime for scales k < kr,. However, 
similar to the renormalized theory, the poles approach 
the VEV limfc_j.o Ps{k) — > Po m the IR limit from below 
and the domain, where W is multi- valued collapses to a 
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FIG. 12: Effective field theory with finite A: Graphical representation of the dimensionless superpotential derivative u(X) as a 
function of X — p — po/k in (45 1 for weak, intermediate and strong superfield coupling (from left to right). Each panel shows 
the RG flow, starting with ua[p) = T {p — po(A)) according to |l9| in the UV up to the IR limit. 



point. Hence, the effective Potential is well defined and 
unique. 

For stronger couplings c < cp, the effective potential 
is plagued by Landau poles and becomes multi-valued 
even in the IR. This becomes apparent by considering the 
second derivative W" of the superpotential. The latter 
shows a non-analyticity at po exactly in the IR limit with 



lim W"{p) = ^—, 

P->-Po± C ± 7T 



(67) 



where W'(p ±) — > ±0. Apparently, the solution W 
shows a cusp with positive W" for W' — > +0 and negative 
W" for W — > —0 in the vicinity of the node if |c| < cp. 
Since there exists at most one Landau pole with W' L < 
in the IR limit (Fig. 11 ) and since W'(p —> — oo) = — oo, 



it becomes apparent that there has to be a Landau pole 
located in the physical regime for k — > if and only if 
Id < c P . 



C. Discussion 

Now we compare and discuss the phase diagrams ob- 
tained by (a) considering the renormalized theory with 
A ^ oo and (b) looking at the effective theory with A 
finite. 

Firstly, let us compare the phase diagrams Fig. [2] and 
[9] as derived from the gap equations (42) and ( JrJTj ). Ap- 
parently, the gap equations (61) of the effective theory 



contain an additional, cutoff (and rcgulator)-dependent 
contribution H(p) compared to (42). The term H(p) 



thus leads to the following modifications of the phase 



diagram of the renormalized theory: In the symmetric 
phase, it diminishes the parameter-range where we ob- 
serve further masses in addition to m. Besides, we find 
up to five different O(N) symmetric phases in the very 
strong-coupling regime |c| <C 1 and for certain VEV pq. 
In the spontaneously broken regime, the function H(p) 
enlarges the parameter range to infinitely large couplings 
r = 1/c, where we observe a second mass M in addition 
to M p . However, since the masses in the very strong 
coupling regime are quite large, i.e. of the order of the 
cutoff A, we believe them to be regulator-dependent and 
unphysical. 

Next, let us compare the phase diagrams Fig.[7|and[l0| 
as deduced from our RG studies. Here, we claimed solu- 
tions W'(p) to be physically relevant, if there exists no 
Landau pole characterized by an infinitely large fermion- 
boson coupling W" in the physical domain. 

Let us first consider the SYM regime. Here, the narrow 
window between the couplings Cl and cp, where there 
exist two masses m and M vanishes for finite A and the 
effective theory shows only a single mass m. Further- 
more, the strong coupling domain is reduced and becomes 
Po-dependent for A finite. The different structure of the 
SYM regimes become apparent by comparing Fig. [3] with 
Fig. 12 13 In the renormalized theory (Fig. [3]), there ex- 
ists an UV Landau pole in the physical domain for super- 
field couplings stronger than t = cj/ and the potential is 
not even defined for all fields p > 0. In contrast, the effec- 
tive theory always features a well-defined UV limit, given 
by the superpotential W' A — t (p — p (A)) at the UV 
scale k = A. The potential is defined for all fields but may 
show ambiguities for sufficiently strong couplings. This 
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FIG. 13: Effective field theory with finite A: Graphical representation of the dimensionful superpotential derivative W' /A as 
a function of X = (p — po) / A for weak, intermediate and strong superfield coupling (from left to right). Each panel shows the 
RG flow of the superpotential, starting with W^(p) — r (p — po(A)) according to ( |19[ | i n the U V up to the IR limit. Note the 



non-analyticity of the effective superpotential at X = in the IR limit k — > (see Sees. VB 2 and VD 1. Note further that the 
running potential with c = cp (middle panel) shows ambiguities for very small \X\ <C 1 for small scales k which are not visible 
in the figure. 



is illustrated in Fig. 13 where the superpotential W'/A 
is plotted as a function of X — (p — po) /A. In the SYM 
phase, the origin p = corresponds to X = |po|/A > 0. 
Now, let us consider the strongly coupled domain with 
\c\ < cp fixed (Fig.[l3j right panel) and |po| "C 1. Here, a 
IR Landau pole occurs at fcj, > in the physical regime 
and additional masses at the origin appear by approach- 
ing the IR. However, if we choose \po\ large enough, the 
IR Landau pole drifts out of the physical domain and 
the effective potential is unique and well-defined for all 
p > with a single mass m. This upper limit of | po | ( c) 
simply corresponds to the borderline connecting (cp,0) 
and (0,-1.077) in Fig. [lO 



We find identical weak, Landau and strong coupling 
SSB regimes for the renormalized and the effective the- 
ory, see Fig. [7J [lO] The renormalized as well as the ef- 
fective theory exhibit an IR Landau pole for all \c\ < cp 
(Fig. [3] and Fig. [T~3] > . The existence of an IR Landau 
pole within the effective theory is shown as follows, see 
Fig. [T3j right panel: The origin p = corresponds to 
X = —p~o/A < and thus there always emerges an 
IR Landau pole in the physical domain at fcj, > for 
|c| < cp. Independent of the superfield coupling and 
the VEV po > 0, there always exists only a single mass 
M at the origin representing an O(N) symmetric phase. 
Again, the effective potential is always defined for all 
fields, but may show ambiguities, whereas the potential 
is not defined for all fields p > in the strong coupling 
regime |c| < cp in the infinite cutoff-limit A — > 00. 



the renormalized and the effective model. Notice that 
these masses represent O(N) symmetric phases of the 
model, since they emerge from the curvature of the po- 
tential at vanishing field p — 0. The parametrically large 
masses m observed in the spike-like region (see Fig. [9]) 
are not included in Fig. 14, since we believe them to be 
an artifact of the chosen regularization. 

In summary, in the SSB phase, and in the symmet- 
ric phase at weak coupling, the difference between the 
(3 >2 )d=3 theory at finite and infinite UV cutoff is minute, 
resulting in equivalent phase diagrams. In the symmet- 
ric phase for c < cp, the difference is more pronounced: 
At finite UV cutoff the fluctuations of the Goldstone 
modes have less "RG time" available to built-up non- 
analyticities in the effective potential. This leads to a 
shift in the effective boundary between weak and strong 
coupling, allowing for a substantially larger domain of 
a regular O(N) symmetric phase. At strong coupling, 
we also conclude that the absence of an O(N) symmet- 
ric phase at infinite cutoff arises from the theory at finite 
UV scale through an O(N) symmetric phase with anoma- 
lously large mass of the order of the UV scale itself. 



D. Effective potential 



Finally, Fig. 14 compares the different mass scales of 



As already mentioned in Sec. |IIIB the relevant micro- 
scopic coupling k — k ci + po/A determines the macro- 
scopic physics of the model: if k < k ci (po < 0), the ef- 
fective potential preserves global O(N) symmetry. Con- 
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FIG. 14: Green, solid lines: mass scales m, M of the renormalized theory (A — > oo) as functions of the coupling c according to 
the gap equations ( |43[ ) for fixed po = {—0.5, —0.5,0.5} (left, middle, right panel). Blue, dashed lines: Masses m/A, M/A of the 
effective theory (A < oo) as functions of the coupling c according to the gap equations p0|) for fixed pa/ A = {—0.5, —0.5,0.5}. 



trary, if n > k c1 . (p > 0), the symmetry may be sponta- 
neously broken, if the VEV po > is taken. The specific 
UV coupling K cr = 1 marks the phase transition between 
the two regimes. Fig. p~5] shows the flow of the effective 
average potential Vk(p) for different values of k, starting 
in the UV at k = A with 



(68) 



according to (19), up to the IR limit k — > 0. Three 



aspects of the potential need to be discussed further: 

Firstly, there exists a strong coupling domain, where 
the effective potential shows ambiguities within the phys- 
ical domain, both in the infinite cutoff limit (Fig. [3]) and 



in the effective theory limit (Fig. 13 right panel). At 
strong coupling, the effective potential admits no physi- 
cal solution for small fields, except for an unphysical one 
with l/\u'\ <C 1 in the effective theory description. This 
result indicates that a description of the theory in terms 
of an effective superpotential is no longer viable, possibly 
hinting at the formation of bound states with or with- 
out the breaking of supersymmetry Incidentally, for the 
same parameter values the effective potential admits two 
solutions for large fields, except for an unphysical third 
solution one in the effective theory description. The the- 
ory admits two different effective potentials associated 
to the same microscopic parameters, which has been dis- 
cussed in pQ in the context of fixed point solutions. 

Secondly, the effective potential at k — is non- 
analytic at its nontrivial minimum p . Consider therefore 
the second derivative of the superpotential 



W" (p) 



1 



W (W'/k) - W (W'/A) 



(69) 



in the vicinity of po(k), where W'(po(k)) — according 
to (59). By approaching the IR, (69) simplifies to (67) 



Apparently, this non-analyticity does not appear until 
the exact IR limit k = is approached. Contrary, for all 
finite scales k > we find W"(po(k)) = 1/c, which sim- 
ply represents the exactly marginal superfield coupling r. 
Since the radial mass is given by 



M 2 p = V"( ( f>)\ (t>=<t>o = (2pW"(p)r\ poiky 



(70) 



a uniquely defined radial mass only exists for finite scales 
k > and reads 

M p (k) = 2 T p Q {k) = 2r{k + po), po > 0. (71) 

First studies at finite N indicate that the non-analyticity 
of W for k — is solely due to the large- N limit. 

Thirdly, the effective scalar field potential in the SSB 
phase with non-vanishing VEV is not convex, even in 
the IR limit k — > 0. As it has already been mentioned in 
|30) , the supersymmetric analogon of the potential term 
in the classical action is the superpotential W, Con- 
sequently, a flow of the superpotential is derived which 
drives the approach to convexity of the superpotential 
W, but not necessarily of the potential V = pW' 2 . The 
superpotential W is a convex function if and only if the 
first derivative W'(p) represents a monotonically increas- 
ing function of p. According to ( 62 ) , this condition is 



satisfied as long as c > cp, i.e. in the weakly coupled 
domain. This fact supports the conjecture that super- 
symmetry may be broken spontaneously in the strongly 
coupled domain exhibiting Landau poles. 
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FIG. 15: RG flow of the effective average potential 14/A 3 as a function of p/A according to (591 for different values of 
5k = po/A = { — 0.1, 0, +0.1} at weak coupling c = 3.7. If po < 0, the system evolves into an O(N) symmetric phase (left 
panel). Vanishing po corresponds to the phase transition between the O(N) symmetric and the SSB phase and the scale 
invariant solution is approached in the IR limit (middle panel). If po > 0, the macroscopic theory is characterized by a 
non-vanishing VEV po(k — > 0) = po > (right panel). The insets show the potential at small fields approaching the IR limit. 



E. Phase transition & critical exponents 

The supersymmetric O(N) model in d — 3 is an effec- 
tive field theory that features the large-distance proper- 
ties of statistical models near a second order phase tran- 
sition. According to [I], the fixed-point solution charac- 
terizing the phase transition shows Gaussian scaling for 
all finite couplings c, except for |c| = cp,cj. Following 
[15) we can also extract the thermodynamical critical ex- 
ponents. The expectation value of the field (</)) serves as 
order parameter, and in the SSB regime it is related to 
the VEV po via (choose (pi = Sn<fi) 

^ifo = VMkA. (72) 



{4>) = lim ^/2p {k) 

k— 



We may associate the deviation of k from its critical value 
n cr = 1 with the deviation of the temperature T from 
the critical temperature T c according to SkA ~ (T c — T). 
Thus we have 



(po)" with = 



1 



(73) 



Next, consider the critical exponent v describing the 
manner in which the correlation length £ diverges (the 
mass vanishes) by approaching the phase transition. We 
thereby distinguish between 

C" 1 = m ~ {-pof (SYM regime, p Q < 0) 

r 1 = m ~ {p f (SSB regime, p > 0). (74) 

Let's consider first the squared masses corresponding to 



O(N) symmetric ground states as given by (41 ). We are 



interested in how the superpotential W' vanishes at the 
origin when po — ¥ 0. We begin with the parameter range 
c > and c ^ cp. Using (62) and (25) we have 



P- Po 



cW 



n\W'\-^W' 2 
A 



O 



fW' 
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V A 3 



(75) 



for small masses. In the SYM regime, this gives (43) for 



W'/A <C 1, where the second mass in (43) only exists in 
the strong coupling region c < n. Hence, according to 
( 74 ) we have 



u = l. 



(76) 



In the SSB regime, there exists a unique O(N) symmetric 
ground state with mass M given by (43) for all c > tt, 
implying 



v' = 1. 



(77) 



We also observe a spontaneously O(N) broken ground 
state, characterized by its radial mass according to ( 71 1 . 



Since M p ~ po, this also leads to (77) 



Now consider the exponent 5, given by J\p o=0 ~ 4> , 
where J = dV/d(f>. Close to the phase transition, where 
we may assume the cutoff to be much larger than the 
mass scale, i.e. W'/A <C 1, the effective potential reads 



V(p) = j^p(p-Po) (78) 
with A = c + 7rsgn(p — po) and sgn(0) = 0. This leads to 

3 ^8 



4A 2 



with 5 



(79) 
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Finally, we discuss the critical exponent 7 associated with 
the susceptibility x — d<j>/dJ — (d 2 V/d(j) 2 )~ 1 near the 
phase transition, 



x(J)l J=0 ~(-Po) 7 
X(J)\ J=0 ~ (Po) Y 



(SYM phase, p < 0) 

(SSB phase, p > 0). (80) 



Using 1781) and (80) we get 



7 = 7 



(81) 



Note that the results (|76j), (|77|, (|79) and (|81j) are in- 
variant under changing c O — c, see (44). The thermo- 



dynamical scaling exponents derived here can equally be 
obtained from the leading RG exponent together with 
scaling relations by using v = 1/6, where 9 = 1 is the IR 
relevant eigenvalue due to the VEV. The scaling expo- 
nents in the special case where c = ±cp are discussed in 
the following section. 
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VI. SPONTANEOUS BREAKING OF SCALE 
INVARIANCE 



FIG. 16: Phase structure of the scalar O(N) model at infinite 
N including the BMB fixed point, according to [5] (see text). 



In this section we discuss the supersymmetric analo- 
gon of the Bardeen-Moshe-Bander (BMB) phenomenon, 
the spontaneous breaking of scale invariance and the as- 
sociated non-classical scaling. 



A. Bardeen-Moshe-Bander phenomenon 



fixed point. One expects that at finite N the tricritical 
line extends all the way to infinite rj and the BMB point 
disappears (43] . We note that the BMB fixed point is also 
of interest as a fundamental UV fixed point, allowing for 
a non-Gaussian continuum limit for the {<fi 2 ) d=3 theory 
with non-classical scaling. 



We first recall the BMB phenomenon for scalar O(N) 
symmetric theories. Linear O(N) models serve as perfect 
testing ground for studying critical phenomena. For large 
N the solvable spherical model gives a qualitatively ac- 
curate picture of the phase structure of the theory. The 
( < t )2 )d=3 theory exhibits an IR-attractive Wilson-Fisher 
fixed point corresponding to a second order phase transi- 
tion between the O(N) symmetric and the spontaneously 
broken phase [20]. In contrast, the scalar {<fi 2 ) d=3 model 
shows a more complex phase structure [3-5 , 20J . Depend- 
ing on the renormalized couplings p 2 , A and rj of the oper- 
ators <fi 2 , <fi 4 and <fi 6 , one observes a first-order phase tran- 
sition without universal behavior or a second-order phase 
transition with universal behavior. Both regimes are sep- 
arated by a tricritical line t, characteri zed by vanishing 
couplings /i 2 and A as depicted in Fig. 
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first-order transitions continues into the 0(N) symmetric 
phase for couplings with rj > r\ c and ends at a gas-liquid 
transition line I. Scale invariance is an exact symmetry of 
the tricritical theory, but at the end point (0, 0, rj c ), scale 
invariance is spontaneously broken. The free coupling 
77 is dimensionally transmuted to an undetermined mass 
scale m and a massless Goldstone-boson (dilaton) shows 
up. In the large-iV limit this non-trivial and UV-stable 
BMB fixed-point marks the point where the tricritical 
line t and the gas-liquid line I meet. Hence, the tricriti- 
cal line connects the Gaussian fixed point and the BMB 



B. Supersymmetric BMB phenomenon 

In the supersymmetric theory, the BMB phenomenon 
has first been discussed in [5] with variational methods. 
Here, the critical (^ 2 ) 2 [ _ 3 theory with a quartic superfield 
potential corresponds, in the scalar sector, to a critical 
( ( t )2 )d=3 with a sextic potential. The main new addition 
due to supersymmetry is that the scalar quartic and sex- 
tic couplings are no longer independent of each other. 

Using the fully integrated RG flow, the following pic- 
ture for the BMB phenomenon emerges: If we fine-tune 
the classical coupling k = k ci , the solution (22) at the 
origin p = reads 



- 1 = cu + H (uq) 



(82) 



in the IR limit, where uo = u(p = 0). This equation 
simply represents the fixed point solution at vanishing 
field. The 0(N) symmetric ground state is characterized 
by the mass 



M' 



= (W\Q)f = p 2 = (u kf . 



(83) 



where M — — ft, > 0. Evidently, uq has to diverge as 
1/k in order to allow for spontaneous breaking of scale 
invariance with a finite mass scale M in the IR limit 
k — >• 0. Now we find the transcendental equation (82) to 
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have always a single zero mass solution M = 0, except 
for c = cp, where it shows an additional, infinitely large 
solution u — > — oo. Note that this limit emerges from 
u(p) through negative field squared values p — > — , which 
is a consequence of our regularisation. Hence, the specific 
microscopic parameters 



(MAO 



(84) 



lead to a macroscopic theory, where the mass of the O(N) 
bosonic and fermionic quanta is left undetermined. Thus, 
scale invariance is spontaneously broken in accordance 
with [5J [TT] and a mass is generated by dimensional 
transmutation. The coupling parameter r takes the value 
( [84] ) in our conventions, and the associated degree of free- 
dom is 'transmuted' to an arbitrary mass scale M. Spon- 
taneously broken scale invariance leads to the appearance 
of a Goldstone boson (dilaton) which is accompanied by a 
Goldstonc fermion (dilatino) , since supersymmetry is left 
unbroken. Note that these particles are exactly massless, 
since r is not renormalized. 



the symmetric and symmetry broken regimes agree. We 
also stress that the BMB scaling exponent (86) is non- 
classical. Furthermore, it cannot be derived from the 
RG scaling alone, as they are due to non-analyticities in 
the field dependences. As a final comment we note that 
an infinite uq, the fingerprint for spontaneous breaking 
of scale invariance, is stable under alterations of the RG 
scheme. 



VII. RADIAL MODE FLUCTUATIONS 

In this section, we give a first account of the phase 
transition in a theory with finitely rather than infinitely 
many supermultiplets N, focussing on the existence of 
a fixed point, the phase transition, and the fate of the 
supersymmetric BMB phenomenon to leading order in a 
gradient expansion. 



C. BMB scaling exponents 



A. Exact fixed point 



Next we turn to the scaling exponents of the super- 
symmetric BMB fixed point. The critical exponents (74) 
and ( 80 ) become double- valued due to a different scaling 
behavior of the different mass scales m, M near the fixed 
point. These, in turn, originate from the finite and the 
infinite uq solutions detected at \c\ = cp, see Fig. [3} The 
latter is responsible for the special nonanalytic behavior 
of the solution at the BMB fixed point. We first con- 
sider c = cp, and the critical exponents v and v' defined 
in (74). By approaching the phase transition from the 



SYM regime, we find m = —pq/2tt and hence 



v = 1. 



(85) 



In turn, approaching the fixed point from the SSB regime, 
the expression for M in (43) is not applicable since the 
contribution linear in W in (75) vanishes. The sub- 



leading quadratic terms take over and we are lead to 
M 2 = fpo, implying that the supersymmetric BMB ex- 
ponent v' is given by 



y BMB 



(86) 



We now consider c = — cp. By virtue of the symmetry 



(44) we note that the mass scales m -H- M interchange 
their roles under cp O — cp. Consequently, the scal- 



ing exponents ( 85 ) and ( 86 ) also interchange their val- 
ues. Therefore we conclude that the theory at |c| = cp 
displays conventional scaling with ( 85 1 as well as un- 
conventional scaling with (86). The former is a conse- 



quence of the smooth 'non-BMB-type' scaling related to 
finite ito, whereas the latter is the BMB scaling associ- 
ated to infinite u . In either case, and under the above 
identification, we conclude that the scaling indices from 



The main new addition to the supersymmetric RG flow 
at finite N are the fluctuations of the radial mode. They 
imply that the quartic coupling r is no longer an ex- 
actly marginal coupling with an identically vanishing (3- 
function. Instead, the flow of this coupling is governed by 
terms of order 1/N. The absence of an exactly marginal 
coupling implies that the line of fixed points found at in- 
finite N will collapse into a finite, possibly empty set of 
fixed points. Furthermore, the running of the VEV no 
longer factorizes from the other couplings of the theory 
resulting in a more complex structure of the RG flow. 

In order to study the supersymmetric O(N) model at 
finite N we return to the full RG flow fi"5"[ ), which in 
terms of u = w' takes the form 



n in 1 n I 1 _ u2 

d t u = -u + pu-(l--)u (TT -^ 

1 , , 9 „ s l-{u + 2pu'f 
-N {6U+2U P \l + {u + 2pu>W 



(87) 



A global, analytical, solution of the RG flow (87) is 



presently not at hand, and we have to resort to approxi- 
mate solutions instead [25] ■ We start with a polynomial 
approximation to order n for the 'potential' u, writing 



i(p,t) =J2a t (t)(p- Pa (t)Y 



It expresses the potential in terms of (n + 1) couplings 
(poj^i)'" , a, n ) to determine its fixed points. Inserting 
the ansatz (88) into the PDE (87) we find a tower of 
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ordinary, coupled differential equations for the couplings, 

1 

A 



dtPo(t) = -po(t) + 1 - 



4p (*)- 



a*(t) \ (l-(2 /O0 (t)q 1 (t)) 2 ) 
ai (i)y (l + (2 / j (t)a 1 (t))2)a 



d t a„(t) = f n (p (t),ai(t),a 2 (t), ■ ■ ■ ,a n+2 (t)) . 



(89) 



Note that the functions /„ depend on the couplings a n +i 
and a n +2 , because the RHS of ( 87 1 involves up to second 
derivatives of u. The fixed point solution requires the 
flow of all couplings to vanish and hence we set the LHS 
of (89) equal to zero, leading to an algebraic system of 
(n + 1) equations for (n + 3) unknowns. These may be 
solved, tentatively, by setting the last two couplings a n +i 
and a n+ 2 to zero. We find 



Po*(N) = 1 
ai*(A) = 



1 

A 
1 A 



27V- 1 
3 N 2 



(90) 



for the first three couplings. The solution bifurcates into 
two independent fixed points starting with 0,3. Intrigu- 
ingly, the recursive relation leads to an exact analytical 
solution of the full system for all N to arbitrarily high ex- 
pansion order n. The reason for this unlikely outcome is 
that the fixed point ( 90 1 is independent of the boundary 



condition which we have imposed initially on the higher 
order couplings. This follows from noticing that all fixed 



point equations ( 89 ) with n > 2 are of the form 







fn(P0i a l> • • • ) a n+2) 

f n (po,ax, . . . ,a„) 

+ (n + 1) (p - 1 + l/N + d t po) a„+i 
n+1 1 - £ 2 



TV 



[(3 + 2n)a n+1 + 2(n + 2)p a n+2 ] 

Q"n+l ■ 



4p £(n + l) 2 (3 ai +4a 2 p )(e 2 -3) 



(1 + e 2 ) 



2VS 



Here £ = 2a\po, and d t po is given according to (891 



At the fixed point (90) we have — 1, and all terms 
proportional to a n+ \ and a n + 2 vanish. Thus, the fixed 
point equation for every a n (n > 2) is independent of 
a n +i and a n+2 provided the first three couplings have 



the values (90), and we are lead to a closed system of 
(n + 1) equations for (n + 1) couplings allowing for an 
exact solution order by order. 



B. Exact scaling exponents 



The new fixed point ( 90 1 has two branches one of which 



is IR attractive in all couplings except for the running 



VEV which remains an IR relevant operator. The second 
fixed point is UV relevant in all couplings and is not 
pursued any further. The universal scaling exponents 
of the Wilson-Fisher type fixed point can be determined 
analytically. From the eigenvalues of the stability matrix 
Bj = d(dtdi)/daj\* we read off that the lowest coupling 
(a = po) defines an IR unstable direction with a critical 
index 



= 1 . 



(91) 



Note that the leading critical exponent v = 1/0q in (91) 
is super-universal and identical to the result at infinite 
N. The exponent does not receive corrections due to 
the fluctuations of the radial mode and therefore cannot 
be used to distinguish universality classes of different N. 
All other couplings dj, i = 1, 2, 3, . . . define IR attractive 
directions with subleading critical exponents 



1 



*(< + !) 




(92) 



The universal eigenvalues 9i are strictly negative for all 
iV > 1. Furthermore, the Gaussian critical exponents 
#G,i = 1 — i for integer i > of the theory in the large- A" 
limit are recovered from (91 1, (92) in the limit l/N — > 0. 
In particular the formerly exactly marginal 4> 6 coupling 
has now become irrelevant. 



Similarly, the fixed-point values of the couplings ( 90 ) 
converge to the large- A fixed-point values for A — > 00. In 
the presence of the radial fluctuations, the A-dependent 
quartic superfield couplings r*(A) is given by the coef- 
ficient ai * (A) , see ( |90| . Taking the limit of infinite N 
singles out a unique value for the quartic superfield cou- 
pling, 



lim rJN) = - 



(93) 



meaning that the line of non-trivial fixed points 
parametrized by the exactly marginal superfield coupling 
r has shrunk to a single point. Notice also that the 
fixed point value ( 93 ) is different from the supersymmet- 



ric BMB value r = 1/cp in the infinite A limit, see (28). 
This serves as a strong indication for the non-existence 
of a supersymmetric BMB fixed point in the presence of 
the radial fluctuations and A > 1. 



C. Global scaling solution 



The infinite A limit ( 93 ) belongs to the strong coupling 



regime where the fixed point solution for the superpoten- 
tial derivative w* displays two branches, neither of which 
extends towards arbitrarily small fields [1 . The latter, 
signalled through the divergence of du* /dp at some finite 
field value p > 0, is responsible for the occurrence of a 
Landau scale. It remains to be seen whether the fixed 
point at finite A continues to belong to the strongly cou- 
pled regime or not. 
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To answer this question, and to compare the fixed 
points at finite and infinite N, we need to study the fi- 
nite N potentials at small fields numerically. The Taylor 
series ( 88 ) of the scaling solution has a finite radius of 
convergence. Alternatively, one may expand the inverse 
fixed point solution p(u) in powers of u. At infinite N, 
the analytical scaling solution p = 1 + cu* + Hiu^) = 
Si=o biul has a finite radius of convergence r set by the 
gap of the inverse propagator (here: r = 1) [H]. Either 
expansion is limited to a finite range in field space. In 
order to cover the full field space, and to make potential 
non-analyticities of the form u'^(p) — > oo visible, we nu- 
merically integrate the differential equation of the inverse 
function p(w*) instead of M*(p). It reads 



=p — u*p — II — 



1 

N 



2\2 



A^ 



12 



2pp"h 



- (u,p' + 2pf 



(94) 



(p' 2 + Kp' + 2p) 2 ) 2 

subject to suitable boundary conditions. The bound- 
ary conditions p(0) = po* and p'(0) — 2po* correspond 



to a singular point of (94) and cannot be used. In- 
stead, we extract boundary conditions for p(u*),p'(u*) 
for \u*\ = 0.01 -C 1 from the polynomial approximation 
to u*(p) of the order n — 9. The combined use of poly- 
nomial expansions and subsequent numerical integration 
is a well-tested technique in critical scalar theories [2"5] . 



Fig. 17 compares the polynomial approximation of the 
scaling solution with the numerical one for N = 3. The 
graph also contains the analytical solution of the the- 
ory at infinite N. The latter is given by the fixed-point 
equation of (87), where we neglect the contribution of 



the radial mode (the term in the second line) and fix the 
free parameter of the solution to ( 93 1 . We find that the 



large- TV solution approximates the finite- N solution very 
well in the vicinity of the node po* an d above, largely 
independently of the chosen value for N > 1. This is en- 



tirely due to the structure of the fixed point ( 90 1 , where 
2oi*/9o* = I- The numerical solutions illustrate further 
that the fixed point solution at finite N shows a similar 
non-analytic behavior characterized by a diverging mass 
term u'„ — > oo, as it appears in the large- N limit for 



strong quartic superfield coupling (cf. Sec. Ill D ) 



We now discuss the A-dependence of the scaling so- 



lution (90). Fig. 18 shows that the fixed point solution, 



displayed for various integer N > 2, always generates a 
diverging du/dp for some positive field values p — p c {N), 
with < p c {N) < po*(N). The solution u*(p) does not 
exist for small < p < p c , for all N considered. Also, 
we find that it*(p c ) becomes increasingly large in mag- 
nitude with decreasing N . Hence, the main effect of the 
competition between the radial mode and the Goldstone 
mode fluctuations, with decreasing N , is a shift of the 
end point p c (N) and the VEV po*(A) towards smaller- 
values. Continuity in N suggests that this pattern per- 
sists for all N > 1 where po* > 0. 

For the supersymmetric Ising model where N = 1, 
the Goldstone modes are absent and the RG dynamics 




P 

FIG. 17: Fixed point solution u t (p) for N = 3. The figure 
compares the polynomial approximation (blue, dashed line) 
with the non-perturbative integration (green, solid line) and 
a large- iV like solution (black, dotted line). 



is controlled by the fluctuations of the radial mode. In 
the limit N —> 1, (90) predicts a vanishing VEV, po = 



and implies the existence of a supersymmetric Ising fixed 
point valid for all fields, though at the expense of a non- 
analytic behavior of it* (p) at vanishing field. Note that a 
direct study of the N — 1 case using the same RG equa- 
tions [31] has also detected a regular Ising fixed point 
analytic in the fields, whose critical eigenvalue 6q = 3/2 
is different from (91). Furthermore, the diverging of all 
higher order couplings ( 90 ) in the limit N 1 together 



with the continuity of the fixed point in N suggests that 
p c — > and |w*(p c )| — > oo in this limit. This behav- 
ior is intriguing inasmuch as the diverging of it*(p — > 0) 
is the fingerprint for the spontaneous breaking of scale 
invariance. It may thus qualify for a novel supersymmet- 
ric BMB phenomenon which originates from the radial 
mode rather than the Goldstone fluctuations. It would 
seem worth to test this picture directly in the supersym- 
metric Ising model without relying on the limit N — > 1 
adopted here. 



To conclude, the fixed point (90) is of the strongly- 



coupled type for all N > 1 as signalled by the same qual- 
itative behavior seen previously at infinite N [T|. Fur- 
thermore, the fluctuations of the Goldstone modes are 
central for the existence of the endpoint in field space 
p c > of strongly-coupled fixed point solutions. At infi- 
nite N, and as a consequence of p c > 0, the phase dia- 
gram at strong coupling is governed by non-analyticities 
at finite RG scales. Due to p c {N) > for N > 1, the 
same type of non-analyticities with an associated Lan- 
dau scale fci control the phase transition associated with 
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FIG. 18: Fixed point solution u*(p) for various N > 1, show- 
ing N = 2, 3, 10, 20 and 100 from left to right (full lines) in 
comparison with the infinite N result (dashed line). With 
increasing N the solutions con ver ges to the exact infinite N 
result with t(N) approaching (93 1. 



the fixed point (90) at finite N. The above behavior at 
strong coupling is thus generic for supersymmetric (<£> 2 ) 2 
theories with N > 1, and to distinguish from the non- 
analyticities at infinite N responsible for, e.g. the con- 
ventional BMB phenomenon. 



VIII. SUMMARY AND CONCLUSIONS 

Analytical solutions of interacting local quantum field 
theories are benchmarks for a deeper understanding of 
concepts and mechanisms in theoretical physics. In this 
work, we have provided a global rcnormalization group 
study of interacting supersymmetric theories in three eu- 
clidean dimensions, the O(N) symmetric ($ 2 ) 2 Wess- 
Zumino theories, continuing a line of research initiated 
in pp. These theories are the supersymmetric versions 
of O(N) symmetric scalar (</> 2 ) 3 theories, which display 
first- and second-order phase transitions, and the seminal 
Bardeen-Moshe-Bander (BMB) mechanism. 

The main new features due to supersymmctry arise 
through the fluctuations of the Goldstone modes, in par- 
ticular at strong coupling, and their competition with 
the fluctuations of the radial mode. In the limit of in- 
finitely many superfields, the radial mode is absent and 
the theory is solved exactly. The phase diagram is then 
controlled by two free parameters, the exactly marginal 
quartic superfield coupling and the vacuum expectation 
value, which takes the role of an infrared relevant cou- 
pling. Locally, the theory has an interacting fixed point 
for all quartic couplings, yet globally the line of fixed 



points terminates at a critical value. At weak coupling, 
the theory displays a second order phase transition be- 
tween an O(N) symmetric and a symmetry broken phase 
with Gaussian scaling, and global supersymmetry re- 
mains intact. At strong coupling, the global effective 
potential becomes multi- valued in certain regions of field 
space, signalled by divergences in the local fermion-boson 
interactions at a finite Landau scale Hl ■ The appearance 
of the characteristic energy scale Ul resolves the long- 
standing puzzle about peculiar degenerate 0{N) sym- 
metric ground states detected previously [SI [7], show- 
ing that these arise, gradually, from the integrating-out 
of strongly-coupled long wave-length fluctuations. In 
this regime, supersymmetry may be spontaneously bro- 
ken. Furthermore, this pattern is largely insensitive to 
whether an infinite or a finite short-distance cutoff is cho- 
sen, solely inducing a shift in the boundary between the 
weakly and strongly coupled regimes. At finite N, and 
to leading order in a gradient expansion, the additional 
fluctuations of the radial mode lift the degeneracy of the 
quartic superfield coupling and the line of fixed points 
collapses to a finite set. Locally, a new Wilson-Fisher 
type fixed point appears with non-Gaussian exponents 
and super-universal scaling in its infrared relevant cou- 
pling. Globally the fixed point belongs to the strongly 
coupled regime, in complete analogy to the strong cou- 
pling behavior observed at infinite N. In its vicinity, and 
with decreasing N, the admixture of radial fluctuations 
shrinks the domain in field space where a Landau scale 
occurs. The scaling solution extends over all fields as 
soon as the Goldstone fluctuations are absent, though at 
the expense of a square-root type non-analyticity in the 
effective potential at vanishing field. 

The availability of a supersymmetric BMB phe- 
nomenon equally depends on the competition between 
Goldstone modes and the radial mode. At infinite N, 
the Goldstone fluctuations lead to the well-known BMB 
fixed point whose scaling exponent v — 1/2 arises due 
to non-analyticities of the infinite N limit. Supersymme- 
try remains intact, and the spontaneous breaking of scale 
invariance leads to the appearance of an arbitrary mass 
scale together with an exactly massless Goldstone boson 
and fermion. The fixed point disappears in the presence 
of both, radial and Goldstone mode fluctuations. The 
BMB mechansim may re-appear provided the Goldstone 
modes are absent altogether, in which case the sponta- 
neous breaking of scale invariance is driven solely by the 
radial mode. A definite conclusion on this point requires 
more study. 

From a structural point of view, the most distinc- 
tive new feature due to supersymmetry at strong cou- 
pling is the build-up of a multi-valued effective poten- 
tial, accompanied by non-analyticities in the polynomial 
interactions at a Landau scale k^. Here, we have es- 
tablished that this phenomenon arises primarily through 
the fluctuations of the Goldstone modes, irrespective of 
whether there are finitely or infinitely many of them. It 
is worth noting that similar non-analyticities have re- 
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cently been observed in the random-field Ising model, 
where the disorder is implemented with the help of Parisi- 
Sourlas supersymmetry [JD]. In these models, the spon- 
taneous breaking of supersymmetry is directly associated 
to the appearance of cusp- like non-analyticities at a finite 
Larkin scale fc^, analogous to the Landau scale found 
here. Provided this similarity persists on a fundamental 
level, it suggests that supersymmetry may be sponta- 
neously broken in the ($ 2 ) 2 theory at strong coupling. 
Conversely, our findings make it conceivable that the oc- 
currence of a Larkin scale is the signature of a multi- 
valued effective potential in disordered Ising models. 

Finally, we stress that the availability of an analytic 
functional RG for supersymmetry was decisive to achieve 
our results, allowing for a controlled and global interpo- 
lation between the short- and long-distance regimes of 
the theory even at strong coupling. It is a virtue of the 
fully integrated RG flow at all scales that the structure 
of the quantum effective theory has become transparent. 
We expect that the combination of analytical and numer- 
ical tools adopted from [25] will prove equally useful for 
the non-perturbative study of supersymmetry in other 
settings and extensions. 
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Appendix A: Conventions 

Relevant symmetry relations and Fierz identities for 
Majorana spinors are = X^> ^l^X — ~)Cf fl ^ an d 
Okfii = — \{pt<y)^-H- One of the main features of the ac- 
tion is its invariance under supersymmetry transforma- 
tions. The latter are characterized by the supersymmetry 
variations <5 £ <^ ,^ , generated by Af — 1 fermionic generator 
Q. We have 

5 e $ l (x) — ie k Qk® l (x) with 
Q k = -id §h - -fififfin , Q k = -id 9k - . (Al) 



Thus, (Al) leads to the supersymmetry variations 



5<j> 1 = dip* = (F l + i0)e and 5F l = iefrf) 4 (A2) 

of the component fields. The anticommuting sector of 
the superalgebra is given by the anticommutator of two 
supercharges 



{Qk,Qi} = 2i 1 Z l B lt . 



(A3) 



The derivation of the supersymmetric flow equation is 
given in appendix B of pQ. 
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